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Abstract— The Casson fluid model, which is very significant in the biomechanics and polymer processing 
industries, is another term used to describe non-Newtonian fluid behavior. This study of Casson fluid model on 
unsteady MHD Casson fluid flow with Soret-Dufour effect past a vertical plate embedded in porous medium in 
the presence of radiation with heat generation/absorption and viscous dissipation is presented in this research 
article as a numerical investigation of non Newtonian Casson fluid with applied effects. Regulating partial 
differential equations have been used to explain the mathematical model of the flow field. The Crank-Nicolson 
implicit finite difference approach has been used to numerically solve non-dimensionalized flow field 
governing equations. Concentration, temperature, and velocity profile effects of non-dimensional factors have 
been investigated using tables and graphs as aids. Tables have also been used to observe fluctuations in factors 
like skin friction, the Nusselt number, and the Sherwood number in relation to other parameters. 


Keywords— Casson Fluid, Magnetohydrodynamics, Order of chemical reaction, Soret and Dufour effects 


Viscous dissipation. 


I. INTRODUCTION 


Shear stress and shear rate have a non-linear relationship that 
non-Newtonian fluid can decipher. Worldwide, non- 
Newtonian solutions are employed in the pharmaceutical and 
chemical industries as well. Examples include oils, 
deodorizers, chemicals, syrups, thick drinks, cleansers, and 
the production of many different colours. In place of 
ketchup, custard, tooth paste, wheat, paint, blood, and 
shampoo, there are still a variety of polymer liquids and salt 
explanations that aren't Newtonian fluids. Shear and shear 
tension are longitudinally correlated in Newtonian fluid. 
Newtonian fluids such as blood, shampoo, soap, certain oils, 
jellies, paints, and other different production & further 
technical needs cannot be represented. Numerous 
researchers, programmers, and scientists have studied 
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various uses. Non-Newtonian fluids are much more difficult 
to investigate than Newtonian fluids, nevertheless. Direct 
evaluation of non-Newtonian fluid properties using Navier- 
Stokes equations is not possible. There are several dynamic 
behaviors of fluid models that have drawn the attention of 
academics, including power law, Bingham plastic, Brinkman 
type, Oldroyd-B, Maxwell, Walter-B, and Jeffrey. The 
Casson fluid model is a non-Newtonian fluid evolutional 
paradigm. Casson developed the Casson fluid to forecast the 
flow behaviors of the pigment-oil solution. Non-Newtonian 
fluids are widely used in engineering fields and other 
industries, and their applications are rather obvious. 
Numerous research on  Non-Newtonian fluid flow 
applications have been conducted in the presence of various 
effects for many years, which has attracted the attention of 
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many mathematicians. Casson fluid, which exhibits elasticity 
in nature and is one of the non-Newtonian fluid kinds, 
includes things like honey, jelly, tomato sauce, and others. 
Casson fluid, which exhibits elasticity in nature and is one of 
the non-Newtonian fluid kinds, includes things like honey, 
jelly, tomato sauce, and others. Casson fluid can also be used 
to treat human blood. Liaquat Ali Lund et al.[1] is 
investigated the magnetohydrodynamic (MHD) flow of 
Casson nanofluid with thermal radiation over an unsteady 
shrinking surface. Shahanaz Parvin et al.[2] are discussed 
effects of the mixed convection parameter, concentration 
buoyancy ratio parameter, Soret-Dufour parameters, and 
shrinking parameter in MHD Casson fluid flow past 
shrinking sheet. Lahmar et al.[3] studied heat transfer of 
squeezing unsteady nanofluid flow under the effects of an 
inclined magnetic field and variable thermal conductivity. 
Mohamed R.Eid et al. [4] investigated numerically for 
Carreau nanofluid flow over a convectively heated nonlinear 
stretching surface with chemically reactive species. Hammad 
Alotaibi et al.[5] introduced the effect of heat absorption 
(generation) and suction (injection) on 
magnetohydrodynamic (MHD) boundary-layer flow of 
Casson nanofluid (CNF) via a non-linear stretching surface 
with the viscous dissipation in two dimensions. Asogwa and 
Ibe [6] investigated numerical approach of MHD Casson 
fluid flow over a permeable stretching sheet with heat and 
mass transfer taking into cognizance the various parameters 
present. Renu et al.[7] assessed the effect of the inclined 
outer velocity on heat and flow transportation in boundary 
layer Casson fluid over a stretching sheet. Ramudu et al. [8] 
highlighted the impact of magnetohydrodynamic Casson 
fluid flow across a convective surface with cross diffusion, 
chemical reaction, non-linear radiative heat. Recently 
Mahabaleshwar et al.[9] discussed the important roles of 
SWCNTs and MWCNTs under the effect of 
magnetohydrodynamics nanofluids flow past over the 
stretching/shrinking sheet under the repercussions of thermal 
radiation and Newtonian heating. Ram Prakash Sharma et al. 
[10] reports MHD Non-Newtonian Fluid Flow past a 
Stretching Sheet under the Influence of Non-linear Radiation 
and Viscous Dissipation. Naveed Akbar et al.[11] 
investigated Numerical Solution of Casson Fluid Flow under 
Viscous Dissipation and Radiation Phenomenon. Elham 
Alali et al. [12] studied MHD dissipative Casson nanofluid 
liquid film flow due to an unsteady stretching sheet with 
radiation influence and slip velocity phenomenon. T. M. 
Ajayi et al. [13] have studied Viscous Dissipation Effects on 
the Motion of Casson Fluid over an Upper Horizontal 
Thermally Stratified Melting Surface of a Paraboloid of 
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Revolution: Boundary Layer Analysis. N. Pandya and A. K. 
Shukla [14] have analyzed Effects of Thermophoresis, 
Dufour, Hall and Radiation on an Unsteady MHD flow past 
an Inclined Plate with Viscous Dissipation. Mallikarjuna B, 
Ramprasad S and Chakravarthy YSK [15] Multiple slip and 
inspiration effects on hydromagnetic Casson fluid in a 
channel with stretchable walls. Bukhari Z, Ali A, Abbas Z, et 
al.[16] The pulsatile flow of thermally developed non- 
Newtonian Casson fluid in a channel with constricted walls. 
Divya BB, Manjunatha G, Rajashekhar C, et al.[17] Analysis 
of temperature dependent properties of a peristaltic MHD 
flow in a non-uniform channel: a Casson fluid model. 
Ahmad Sheikh N, Ling Chuan Ching D, Abdeljawad T, et al. 
[18] A fractal-fractional model for the MHD flow of Casson 
fluid in a channel. Haroon Ur Rasheed, Saeed Islam, 
Zeeshan, Waris Khan, Jahangir Khan and Tariq Abbas [19], 
Numerical modeling of unsteady MHD flow of Casson fluid 
in a vertical surface with chemical reaction and Hall current. 


Our goal is to shed light on the impact of a vertical porous 
plate with Soret-Dufour, radiation, heat 
generation/absorption source/sink, and higher-order chemical 
reaction in this inquiry. With the aid of tables and figures, 
the impact of different physical parameters on velocity, 
temperature, and concentration profiles is explained. On the 
other hand, tables are used to discuss crucial physical 
parameters like shearing stress, the Nusselt number, and the 
Sherwood number. 


П. MATHEMATICAL MODELING 


The unsteady MHD viscous, incompressible electrically 
conducting fluid's casson flow past an impulsively begun, 
infinitely inclined, porous plate with changeable temperature 
and mass dispersion has been taken into consideration. The 
plate is placed in a porous material and is vertical. The x-axis 
is considered perpendicular to the plate, and the y-axis 
parallel to it. Additionally, it is first believed that the 
radiation heat flux in the x-direction is much smaller than 
that in the y-direction. The fluid's temperature and 
concentration are the same for the plate. The plate's 
temperature and concentration fall exponentially as a result 
of the impulsive motion along the x-axis against the 
gravitational field with constant velocity uo at time t. Since 
the induced magnetic field is very small and the transversely 
applied magnetic field's magnetic Reynolds number is also 
very small, it can be regarded as inconsequential. Cowling 
[21], the flow variables are just functions of y and t since the 
x-direction is infinite. For this problem with an unstable flow 
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field, the governing partial differential equations are given by: 


Continuity equation: 


eus di => y = —v (constant ) (1) 
ду 

2.1 Momentum equation: 

= ЕЕ pam _ Е _ _ B^ ‚2 = 
СР Ре 9 9 LEET-T,)e лде e Lcd Q) 
Ot ду poy (19m ) pK 


2.2 Energy equation: 


Г ої T D,K, ôC Y (mox 
ps C, CPT А EQ e m re yu 2 -0,(T -T,) P 
Ot дуў ду ду c, ду y 


BC OC OC D,K, OT P 
Cy hop mK, OT ктк à) 
Е = —2 —2 r оо 
дї ду ду T, Oy 
Fan 1. n . + + . + 
k C-C terms in mass equation for higher order chemical reaction 
r о q 8 
n order of chemical reaction 
k ; chemical reaction constant 
C concentration 
T temperature 
Т. temperature of free stream 
C, concentration of free stream 
В Casson parameter 
Вс coefficient of volume expansion for mass transfer 
Br volumetric coefficient of thermal expansion 
Tm mean fluid temperature 
q, radiative heat along y * -axis 
Q, Coefficient of heat source/sink 
v kinematic viscosity 
K coefficient of permeability of porous medium 


Dm molecular diffusivity 
k 


thermal conductivity of fluid 
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Cp specific heat at constant pressure 
u viscosity 

p fluid density 

o electrical conductivity 

g acceleration due to gravity 

Кт thermal diffusion ratio 

m Hall current parameter 


NÉ S n 
In Equation(4); k " (C -C m ) has come on account of п? order chemical reaction. 


The boundary conditions for this model are assumed as: 


ї<0; и=0, еу, Pel, С=С, V y 
i50; u=u, T=T,+(T,-T, Je“, C-C,4(C,-C,)e" at y=0 (5) 
u — 0, T 9T, C 2C, as yo 
2 
Where A — 9 
v 


Roseland explained the term radiative heat flux approximately as 
T4 
4o, oT 
B = 
3a,, OY 


q, = (6) 


Here Stefan Boltzmann constant and absorption coefficient are О’, and Я, respectively. 


' ' d T 4 : : 
In this case temperature differences are very-very small within flow, such that T ' canbe expressed linearly with temperature. It 
is realized by expanding in a Taylor series about Т. and neglecting higher order terms, so 


T 247 T-S3T- (7) 
With the help of equations (6) and (7), we write the equation (3) in this way 
T T T 167, UD D,K, ӘС 
Bac, OF yt m d 671. с, © n a m ОС 
Ot Oy 


=) L2 2 
Оу 3a, Oy B. Oy (8) 


+22) -бт-т.) 


ду 


Let us introduce the following dimensionless quantities 
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px usn. ж = quede ge Д 
и, V V T,-T, C cc. 
BD TT. P 
G, = SPIE Co) g 2" 85 E 2) K=% K, 
Uo Vo Uo Vo V 
=a БО 
5. =, poft В = E puis p a (9) 
D k kk Tv(C.—C.J 
D K.(C.-—C,. : 
Q--&V. p, — "m (С, Co) E ДЕЕ 
ОС , Vo C,C,V( T, =]; ) Vo » 
2 " 
a _ oB, v sin’ т 


Vo” " 
== == Й i == 2 5 m 
c,( T, - T, ) (1+ т? jv, 


Using substitutions of Equation 9, we get non-dimensional form of partial differential Equations 2, 8 and 4 respectively 


2 
с jx 2 £+6,046,C-[M, e e (10) 
Ot Oy p joy K 
2 2 : 

00 00 1 (E)n, Sns, CLAR” "m 
дї ду Р 3 ) dy ду ду 

2 2 
C C 1C .00 „ү T 
Ot Oy S, ду” ду” 


With initial and boundary conditions 


ї<0; и=0, 0-0, С=0 Vy 
ї> 0; и=1, O=e", C=e' at y=0 
и — 0, Ө 0, C0 as y о 


(13) 


The degree of practical interest includes the Skin friction coefficients т, local Nusselt Nu, and local Sherwood Sh numbers are 


given as follows: 
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1 | ðu 00 
т=—|1+— == 


ди [oc 
8 ду у= и 


, Sh = –| — 
ду у=0 


— (14) 
OY)y =0 


III. NUMERICAL METHOD OF SOLUTION 


Exact solution of system of partial differential Equations 10, 11 and 12 with boundary conditions given by Equation 13 are 
impossible. So, these equations we have solved by Crank-Nicolson implicit finite difference method. The Crank-Nicolson finite 
difference implicit method is a second order method in time (o(At?)) and space, hence no restriction on space and time steps, that 
is, the method is unconditionally stable. The computation is executed for Ay = 0.1, At = 0.001 and procedure is repeated till y = 4. 
Equations 10, 11 and 12 are expressed as 


- Е -2u. .+ + -2 t 
"pb t PIT H ызгы uu ЖЫ кы unb eb dud 


At Ay 8 2( Ay )’ (15) 
8. 478 C; te DL ы! | 
2 2 К 2 
б жү dg М. d S qi ch teu DUM Lg а у id 
At Ay PA 3 2( Ay f 
ЕК ш ove quoto Жош ш ПОТ. ТЫТ 
+р, 2 u^ 
2( Ay ) Ay 2 
(16) 
a aj ud К шш eg eee ш und 
At Ay Sc X Ay) 17) 
dg urea aut м мм | S pt 
t$, ; tK, 
X Ay) 2 
Initial and boundary conditions are also rewritten as: 
а ЕО: . n =0, = Vi 
"i0 О 9; 0 О Cro O i 
Е 7 | (18) 
. = 1, Ө . = au C . = pe V 
"o, у O, j e 0, j e J 
и, .—0, 0,.50, C, . 0 
Lo l, j L, j 
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Where index 1 represents to y and j represents to time t, At — 
tj*1-tj and Ay = yi+l-yi. Getting the values of u, 0 and C 
at time t, we may compute the values at time t+At as 
following method: we substitute i = 1, 2, ...,  -1 , where n 
correspond to oo ,equations 15 to 17 give tridiagonal system 
of equations with boundary conditions in equation 18, are 
solved employing Thomas algorithm as discussed in 
Carnahan et al.[20], we find values of 0 and C for all values 
of y at t + At. Equation 15 is solved by same to substitute 
these values of 0 and C, we get solution for u till desired 
time t. 


IV. ANALYSIS ON OBTAINED RESULTS 


The present work analyzes the boundary layer unsteady 
MHD Casson flow past a porous vertical plate with the 
Soret-Dufour effect and Hall current. The influence of the 
order chemical reaction has been incorporated in the mass 
equation. In order to see a physical view of work, numerical 
results of velocity profile u, temperature profile 0, 
concentration profile C have been discussed with the help of 
graphs and skin friction coefficients, Nusselt number and 
Sherwood number are discussed with the help of tables. The 
following values are used for investigation Gr = 4.2, Gm = 
6, К = L5, Mi- 02, B= 0.35, Kr = 14, Pr = 0.6, Du = 0.2, 
Sc = 0.25, Sr =1.7, К =1.8, Ec =2, Q=3,t=0.1. 


It is noted from figure 7 that increasing radiation parameter 
R, velocity u increases. This is correct observation because 
the increase in radiation reveals heat energy to flow. It is 
analyzed that an increase in R, temperature 0 increases and 
it is notable that an increase in R, concentration C near to 
plate decrease after that increases in figure 25. In figure 6, 
velocity decreases as Prandtl number Pr increases and 
temperature decreases in figure 15 when Pr increases. In 
figure 21 concentration C near to plate increases and some 
distance from plate concentration decreases as Prandtl 
number increases. Figure 16 depicts the importance of 
radiation on temperature distribution. Figure 23, depicts the 
variation of Schmidt number Sc as concentration decreases 
rapidly with increase Sc while velocity profile in figure 9 
decreases near to plate. In figure 1, 19 and 12, it is seen that 
velocity increases and concentration decreases as increase 
Dufour number Du, whereas temperature increases as Du 
increases. Figure 3 shows that increment in porosity 
parameter K results in an increase in velocity. Figures 8 and 
22 depict the behavior of chemical reaction parameter Kr on 
velocity and concentration respectively. It is seen that 
velocity decreases, concentration decreases rapidly as Kr 
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increase. The negative value of Q < 0 means heat 
absorption and the positive value of О > 0 means heat 
transfer. In figure 4, velocity profile decreases on increasing 
heat source/sink parameter Q and also reducing momentum 
boundary layer. Figure 13 analyzed the impact of heat 
source/sink parameter Q in the temperature profile. It can be 
seen that the temperature profile decreases rapidly and the 
thermal boundary layer reduces for an increase of heat 
source parameter but it increases with the heat sink 
parameter. Figure 20, depicts that concentration profile 
increases near to plate from middle of boundary layer it 
decreases as well as species boundary layer reduces on an 
increase of heat source/sink parameter. Figure 11, 18 and 26 
that velocity, 
increase on increase of time. Figure 14 described that 


reveals temperature and concentration 
increment in Ec results in temperature increases. In Figure 
2, velocity near plate increases rapidly and then decreases 
rapidly. Figure 5 reveals that an increase in hall current 
parameter Му then velocity increases slowly. Figure 10, 24 
and 17 depict that increment in Soret number Sr, velocity 
and concentration increases while temperature decreases 
respectively. 


It is observed from Table 1 that Change in Schmidt number 
Sc effects as skin friction coefficient and Sherwood number 
increases while Nusselt number decreases. Pr effects as skin 
friction coefficient and Sherwood number decrease while 
Nusselt number decreases. Skin friction coefficient and 
Sherwood number Sh increase whereas Nusselt number 
decreases with Dufour number Du increases and Eckert 
number Ec increases. Increase in Soret number Sr, Skin 
friction and Nusselt number Nu increase while Sherwood 
number Sh decreases. On increasing Casson fluid parameter 
p results in skin friction coefficient and Sherwood number 
Sh increases and Nusselt number decreases. On increasing 
Hall current parameter M; results in skin friction coefficient 
and Sherwood number Sh increases and Nusselt number 
decreases. It is also noted that increment in Q heat 
source/sink the skin friction coefficient and Sherwood 
number decrease while Nusselt number increases. 


V. CONCLUSION 


Effect of Hall current, viscous dissipation, first order 
chemical reaction, changes in Soret-Dufour effects on 
unsteady MHD flow past a vertical porous plate immersed 
in a porous medium are analyzed. This investigation the 
following conclusions have come: 


Page | 11 


Shukla 
International Journal of Chemistry, Mathematics and Physics (UCMP), Vol-6, Issue-4 (2022) 


5.1 It has been observed that on increasing Hall current 121 
parameter the velocity of the fluid increases. | 
5.2 Тһе effect of radiation on concentration is noteworthy. 1.0% K = 0.5, 2.3, 4.2, 6.3 
It is observed that increasing values of R, concentration asl 
falls down and after some distance from the plate, it [ 
goes up slowly-slowly while velocity and temperature "asl 
increases. [ 
5.3 Interestingly, increment in concentration has been sa 
found on increasing Prandtl number P r. 02! 
5.4 For increasing values of Kr, it is a considerable | 
enhancement іп velocity, i.e. velocity decreases slowly 0.07 DENEN NEAN o "M 4 


but concentration decreases rapidly. 

5.5 Increasing values of Dufour number, it is observed that 
velocity and temperature profile in the thermal Fig. 3: Velocity Profiles for Different Valuesof K 
boundary layer increases whereas concentration profile 
first decreases after then increases slowly in the [ ' ' — 
boundary layer. L 

5.6 Schmidt number greatly influences the concentration | 
profile in the concentration boundary layer. 


12 F 


10} Du = 0.2, 3.2, 5.6, 8.2 
0.8} 
= 0.6} 
0.4} 
0.2} 
0.0L 
0 1 2 3 4 
y 12[ 
Fig. 1: Velocity Profiles for Different Values of Du [ 
1.0% М, = 0.2, 2.1, 4.5, 6.3 
r ] 0.8 
1.2 ] | 
[ ] а 
f 0.6} 
ка В = 0.2 2.3, 5.2,8 | 
[ 0.4 
0.8 [ 
"E 02[ 
0.6} | 
0.0. : 
0.4 0 1 2 3 4 
L y 
0.2 
L Fig. 5: Velocity Profiles for Different Values of M; 
0.0L 


y 


Fig. 2: Velocity Profiles for Different Values of 8 
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1.2} | 


Pr= 0.5, 3.2, 5.5, 3.9 


1.0} 
0.8} 
0.6 | 


0.41 


0.2} 


0.0L 
0 


Fig. 6: Velocity Profiles for Different Values of Pr 


К = 0.2, 2.8, 4, 7 


Fig. 7: Velocity Profiles for Different Values of R 


Г Kr = 1.5, 4, 7.2, 9.5 


Fig. 8: Velocity Profiles for Different Values of Kr 
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1.2} 
Sc = 0.25, 1.75,2.8, 4.2 


1.0} 
0.8} 
0.6} 
0.4 


02. 


0.0L 
0 


10 Sr = 0.2, 1.7, 2.5, 4.7 


Fig. 10: Velocity Profiles for Different Values of Sr 


t=0.1, 0.2, 0.3, 0.4 


“0 1 2 3 
т 
Fig. 11: Velocity Profiles for Different Values of t 
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[ " r T r r r 7——T " ы т 


Du = 0.2, 3.2, 5.6, 82 


9=—3.2, -0.3,1, 3, 5, 7 


1.2 Ес = 1, 2.5, 4, 7 


Fig. 14: Temperature Profiles for Different Values of Ec 
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0.8 


Pr = 0.5, 3.2, 5.5, 8.9 


0.6 


0.4 
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Fig. 15: Temperature Profiles for Different Values of Pr 


Г” " r " T r r r 7—1 " ы " ы т 


R-0.2,2.8,4,7 


"X Кылыз сы лы E ыы RET MEI 


Sr = 0.2, 1.7, 2.5, 4.7 


Fig. 17: Temperature Profiles for Different Values of Sr 
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t=0.1, 0.2, 0.3. 0.4 


1.0} Pr = 0.5, 3.2, 5.5, 8.9 


“0 1 2 3 4 0 
Y y 
Fig. 16: Temperature Profiles for Different Values of t Fig. 21: Concentration Profiles for Different Values of Pr 


Kr — 1.5, 4, 7.2, 9.5 


12 ] [ 
Du = 0.2, 3.2, 5.6, 8.2 ] 0.8 - 


y 


Fig. 22: Concentration Profiles for Different Values of Kr 


12 ] 
Sc = 0.25, 1.75,2.8, 4.2 


ior \ Q-—32, -0.3, 1, 3, 5, 7 | is 


Fig. 23: Concentration Profiles for Different Values of Sc 


Fig. 20: Concentration Profiles for Different Values of Q 
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T T T T ы r T T T T T T T T T T T T T T T 


12 


Lof N Sr = 0.2, 1.7, 2.5, 4.7 


t=0.1, 0.2, 0.3, 0.4 


0.6} 
0.4} 


0.2L 


0.0L 
0 


Fig. 24: Concentration Profiles for Different Fig. 26:Concentration Profiles for Different Values of t 
Values of Sr 


К = 0.2, 2.8, 4,7 


0.0L imi 
0 1 2 3 4 
y 


Fig. 25: Concentration Profiles for Different Values of R 


Table 1. Skin friction coefficient t, Nusselt number Nu and Sherwood number Sh for different values of parameters taking fix 
values of Gr = 4.2, Gm = б, K = 1.5, Mi- 0.2, B= 0.2, Kr = 1.4, Pr = 0.6, Du = 0.2, Sc = 0.25, Sr =1.7, R =1.8, Ec =2, О =3, t 


= 0.1. 
В Du | Ec |K |Kr Pr R Mi |$ Sr Q t T Nu Sh 
02 |02 |2 15 |14 | 0.6 1.8 [0.2 | 0.25 | 17 3 | 0.1 1.47062 1.0054 0.873165 
23 | 0.2 | 2 15 |14 | 0.6 1.8 [0.2 | 0.25 | 17 3 | 0.1 1.55363 0.948078 | 0.889804 
52 |02 |2 15 |14 | 0.6 1.8 [0.2 |025 | 17 3 |01 1.48919 0.937503 | 0.892974 
8 |02 |2 15 |14 | 0.6 18 [0.2 | 0.25 | 17 3 |01 1.46724 0.93406 0.894012 
0.35 | 0.2 | 2 15 |14 | 0.6 1.8 [0.2 |025 | 17 3 |01 1.65241 0.992201 | 0.876917 
035 | 3.2 |2 15 |14 | 0.6 1.8 [0.2 Á |025 | 17 3 |01 1.85239 0.764003 | 0.947336 
0.35 | 5.6 | 2 15 |14 | 0.6 1.8 [0.2 |025 | 17 3 [0.1 |2.03522 0.514929 | 1.02852 
0.35 | 82 | 2 15 |14 | 0.6 1.8 [0.2 | 0.25 | 17 3 |0.1 |2.26875 0.109025 | 1.16889 
0.35 |02 1 |15 |14 | 06 1.8 [0.2 | 0.25 | 17 3 [01 1.62226 1.02237 0.868285 
0.35 |02 | 25 | 15 | 14 | 06 1.8 [0.2 | 0.25 | 17 3 |01 1.66744 0.977248 | 0.881188 
0.35 |02 4 |15 |14 |06 1.8 |0.2 |025 | 17 3 |01 1.71233 0.932884 | 0.89383 
0.35 |02 7 |15 |14 |06 1.8 [0.2 | 0.25 | 17 3 |01 1.80132 0.846208 | 0.918412 
0.35 |02 |2 05 | 14 | 0.6 1.8 [0.2 | 0.25 | 1.7 3 | 0.1 | 0.950611 0.992346 | 0.877565 
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035 |02 |2 2.3 |14 | 0.6 1.8 |02 | 0.25 1.7 3 0.1 1.78111 0.991901 0.876891 
0.35 | 0.2 | 2 42 |14 | 0.6 1.8 | 0.2 | 0.25 1.7 3 0.1 1.89192 0.991573 0.876892 
0.35 | 0.2 | 2 6.3 | 1.4 | 0.6 1.8 | 0.2 | 0.25 1.7 3 0.1 1.93707 0.99142 0.876899 
0.35 | 0.2 | 2 1.5 1.8 | 0.6 1.8 | 0.2 | 0.25 1.7 3 0.1 1.64311 0.991866 0.889224 
0.35 | 0.2 | 2 1.5 4 0.6 1.8 | 0.2 | 0.25 1.7 3 0.1 1.43161 0.983849 1.17117 
0.35 | 0.2 | 2 1.5 7.2 | 0.6 1.8 | 0.2 | 0.25 1.7 3 0.1 1.21004 0.974585 1.47333 
0.35 | 0.2 | 2 1.5 9.5 | 0.6 1.8 | 0.2 | 0.25 1.7 3 0.1 1.07724 0.968551 1.65972 
0.35 | 0.2 | 2 1.5 1.4 0.5 1.8 |02 | 0.25 1.7 3 0.1 1.71762 0.911549 0.901144 
0.35 | 0.2 | 2 1.5 1.4 3.2 18 |02 | 0.25 1.7 3 0.1 1.1256 2.15668 0.470401 
035 |02 |2 1.5 14 5.5 1.8 |02 | 0.25 1.7 3 0.1 0.993012 2.78179 0.229897 
0.35 | 0.2 | 2 1.5 1.4 8.9 1.8 |02 | 0.25 1.7 3 0.1 0.891104 3.51046 -0.059222 
0.35 | 0.2 | 2 1.5 14 | 0.6 0.2 | 0.2 | 0.25 1.7 3 0.1 1.32212 1.56809 0.685784 
0.35 | 0.2 | 2 1.5 14 | 0.6 2.8 | 0.2 | 0.25 1.7 3 0.1 1.77127 0.850637 0.918876 
0.35 | 0.2 | 2 1.5 1.4 | 0.6 4 0.2 | 0.25 1.7 3 0.1 1.87712 0.742687 0.948992 
0.35 | 0.2 | 2 1.5 14 | 0.6 7 0.2 | 0.25 1.7 3 0.1 2.05248 0.593058 0.988149 
035 |02 |2 1.5 1.4 | 0.6 1.8 0.2 | 025 1.7 3 0.1 1.65241 0.992201 0.876917 
035 |02 |2 1.5 1.4 | 0.6 1.8 2.1 | 025 1.7 3 0.1 2.51597 0.988477 0.87732 
035 |02 |2 1.5 14 | 0.6 1.8 4.5 | 0.25 1.7 3 0.1 2.67912 0.987309 0.877551 
0.35 | 0.2 | 2 1.5 1.4 | 0.6 1.8 6.3 | 0.25 1.7 3 0.1 2.70604 0.987102 0.877594 
035 |02 |2 1.5 14 | 0.6 1.8 | 0.2 0.25 1.7 3 0.1 1.65241 0.992201 0.876917 
0.35 | 0.2 | 2 1.5 14 | 0.6 1.8 | 0.2 1.75 1.7 3 0.1 0.579633 0.95069 2.06938 
0.35 | 0.2 | 2 1.5 14 | 0.6 1.8 | 0.2 2.8 1.7 3 0.1 0.360793 0.933651 2.58939 
035 |02 |2 1.5 1.4 | 0.6 18 |02 4.2 1.7 3 0.1 0.18523 0.913773 3.18956 
035 |02 |2 1.5 1.4 | 0.6 1.8 |02 | 0.25 0.2 3 0.1 1.43314 0.985206 1.04842 
0.35 | 0.2 | 2 1.5 14 | 0.6 1.8 | 0.2 | 0.25 1.7 3 0.1 1.65241 0.992201 0.876917 
0.35 | 0.2 | 2 1.5 14 | 0.6 1.8 | 0.2 | 0.25 2.5 3 0.1 1.77073 0.995712 0.783784 
0.35 | 0.2 | 2 1.5 14 | 0.6 1.8 | 0.2 | 0.25 4.7 3 0.1 2.10084 1.00453 0.522318 
0.35 | 0.2 | 2 1.5 14 | 0.6 1.8 | 0.2 | 0.25 1.7 -3.2 | 0.1 -2.01681 0.251471 1.07677 
0.35 | 0.2 | 2 1.5 14 | 0.6 1.8 | 0.2 | 0.25 1.7 -0.3 | 0.1 - 1.82303 0.638556 0.976345 
0.35 | 0.2 | 2 1.5 14 | 0.6 1.8 | 0.2 | 0.25 1.7 1 0.1 -1.7504 0.787159 0.93551 
0.35 | 0.2 | 2 1.5 14 | 0.6 1.8 | 0.2 | 0.25 1.7 3 0.1 -1.65241 0.992201 0.876917 
0.35 | 0.2 | 2 1.5 14 | 0.6 1.8 | 0.2 | 0.25 1.7 5 0.1 - 1.56829 1.17375 0.822771 
0.35 | 0.2 | 2 1.5 14 | 0.6 1.8 | 0.2 | 0.25 1.7 7 0.1 - 1.49549 1.33633 0.772438 
0.35 | 0.2 | 2 1.5 14 | 0.6 1.8 | 0.2 | 0.25 1.7 3 0.1 | -0.759763 0.993876 1.00151 
035 |02 |2 1.5 1.4 | 0.6 1.8 |02 | 0.25 1.7 3 0.2 1.65241 0.992201 0.876917 
0.35 | 0.2 | 2 1.5 14 | 0.6 1.8 | 0.2 | 0.25 1.7 3 0.3 3.25035 1.0526 0.883938 
0.35 | 0.2 | 2 1.5 14 | 0.6 1.8 | 0.2 | 0.25 1.7 3 0.4 4.59795 1.12637 0.938618 
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